function y = e_beale(x)
%  Extended Beale Function  U63 (MatrixRom)
%x0=[1, 0.8, ...., 1, 0.8]
%fmin=n/2
%Author:shi yangyang
n=length(x);
m=n/2;
y=ones(1,m);
for i=1:1:m
    t1=1.5-x(2*i-1)+x(2*i-1)*x(2*i);
    t2=2.25-x(2*i-1)+x(2*i-1)*x(2*i)*x(2*i);
    t3=2.625-x(2*i-1)+x(2*i-1)*x(2*i)*x(2*i)*x(2*i);
    y(i)=t1*t1+t2*t2+t3*t3;
end
y=sum(y);